Distributions on unbounded moment spaces and random 

moment sequences 



Jan Nagel Holger Dette 

^^ \ Ruhr-Universitat Bochum Ruhr-Universitat Bochum 

(N ■ 

Fakultat fur Mathematik Fakultat fur Mathematik 

44780 Bochum, Germany 44780 Bochum, Germany 

e-mail: jan.nagel@rub.de e-mail: holger.dette@rub.de 

p^ 1 July 21, 2010 

Oh' 

Abstract 

In this paper we define distributions on moment spaces corresponding to measures on 
the real line with an unbounded support. We identify these distributions as limiting distri- 
butions of random moment vectors defined on compact moment spaces and as distributions 
corresponding to random spectral measures associated with the Jacobi, Laguerre and tier- 



ce 



mite ensemble from random matrix theory. For random vectors on the unbounded moment 
spaces we prove a central limit theorem where the centering vectors correspond to the 



o 
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moments of the Marchenko-Pastur distribution and Wigner's semi-circle law. 
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1 Introduction 

For a set T C K. let V(T) denote the set of all probability measures on the Borel field of T with 
existing moments. For a measure \x G V(T) we denote by 

mfc(yu) = / x k fi(dx) ; k — 0, 1, 2, . . . 



the k-th moment and define 
(1.1) M(T) = {m{ft) 



(mi(/i),m 2 (/x), 



fieV(T)}c 



as the set of all moment sequences. We denote by II n (n G N) the canonical projection onto the 
first n coordinates and call 



;i.2) 



M n {T) = Ii n {M{T))(Z 



the n-th moment space. Moment spaces of the form (11.11) and (11.21) have found considerable inter- 
est in the literature [see lKarlin and Studdenl ( 119661 )]. Most authors concentrate on the "cl a ssical " 
moment space correspondin g to measures on t he interval [0 , 1] [see Karlin and Shapeleyl fjl953h , 



Krein and Nudelman 



(J1977I ). among others]. IChang et al.l ( 119931 ) equipped the n-th moment 



space M. n ([0, 1]) with a uniform distribution in order to understand more fully its shape and 
the structure. In particular, these authors proved asymptotic normality of an appropriately 
standardized version of a projection Uk(m n ) of a uniformly distributed vector rrin on .M n ([0, 1]). 
Since this seminal paper, sev e ral au thors have extended these investigations in various directions. 
Gamboa and Lozada-Chand ( 2004 ) con sidered large de v iation principles for random moment se- 
quences on the space .M n ([0, 1]), while iLozada-Changl ( 120051 ) investigated similar problems for 
moment spaces correspondin g to more general functions defined on a bounded set. More recently, 
amboa and Rouaultl (J2009l ) discussed random spectral measures related to moment spaces of 
measures on the interval [0, 1] and moment spaces related to measures defined on the unit circle. 

The present paper is devoted to the problem of defining probability distributions on unbounded 
moment spaces. We will investigate these distributions from several perspectives. In Section 2 
we introduce a class of general distributions on the moment space corresponding to measures 
defined on a compact interval. By a limiting argument we will derive canonical distributions 
on the moment spaces corresponding to measures on the unbounded intervals [0, oo) and R, 
respectively. In Section 3 we show that these distributions appear naturally in the study of 
random spectral measures of the classical Jacobi, Laguerre and Gaussian ensemble. Finally, in 
Section 4 we consider random moment sequences distributed according to the new probability 
distributions on the unbounded moment spaces. In particular, we prove weak convergence of a 
centered random moment vector, where the centering vector corresponds to the moments of the 
Marchenko-Pastur law (in the case of the moment space .A/f ([0, oo))) and to the semi-circle law 
(for the moment space Ai(M.)). 



2 Distributions on unbounded moment spaces 

2.1 Canonical moments and recurrence coefficients of orthogonal poly- 
nomials 



Chang et al.l ( 119931 ) considered random vectors on the moment space .M„([0, 1]) governed by a 



uniform distribution. In the present section we will define a class of more general distributions on 
the n-th moment space A4 n ([a,b]) corresponding to the set V([a,b]) of all probability measures 
on the interval [a, b]. The motivation for considering this class is twofold. One the one hand 
we want introduce distributions on the moment space Ai k ([a,b]), which are different from the 
uniform distribution. On the other hand we want to define distributions on unbounded moment 
spaces as limits of distributions on Ai k ([a,b]), when b — a — > oo. For these purposes we will 
make extensive use of the theory of canonical moments, which is briefly defin ed here for the sake 
of a s elf contained presentation. For details we refer to the monograph of IDette and Studden 
( 1199 70 . Let m k -x =(mi, . . . ,m k _i) T E Ai k -i([a, b}) be a given vector of moments of a probability 
measure on the interval [a, b], then these first k — 1 moments impose bounds on the k-th moment 
m k such that the moment vector m k = (mi, . . . ,m k _i,m k ) T is an element of the k-th moment 
space Ai k ([a,b}). More precisely, define for rrik-i E A / ffe_i([a, b]) 

b 

fdfi(t) — rrii for i — 1, . . . , k — 1 



m. 



m. 



mm 



max 



m k (n) 



m k (/j) 



fieV([a,b]) with 
H EV([a,b]) with 



tfd{i(t) = vfii for i = 1, . . . , k — 1 



then it follows that m k = (mi, . . . , m k 
Int C denotes the interior of a set Cstf 
the canonical moment of order / = !,... 



T E Int M.k{[ a i b}) if and only if m k < m^ < m^, where 
Consequently, we define for a point ra k E Int J\Ak([a, b]) 



pi =pi{m k ) 



mi — mf 



I 



k. 



m, 



(0 



,i); l 



k; and that p k describes 



(2-1) 

" l l ~ "H 

Note that for ra^ E Int Aik([a, b]) we have pi E 

the relative position of the moment m^ in the set of all possible k-th moments with fixed 

moments m 1; . . . ,771^-1- It can also be shown that the canonical moments do not depend on 

the interval [a, 61, that is they are invariant under linear transformations of the measure [see 

Dette and Studden! (J1997I )]. Moreover, the definition (12.11) defines a one-to one mapping 

Int Mt b] — >(0,l) n 
m n ^ p n = (p 1 ,...,p n ) T 

from the interior of the moment space Ain ' onto the open cube (0, l) n . It can be shown that 
for a point (mi, . . . ,m2 n -i) E Int A^2n-i([a, b]) the canonical moments appear in the three-term 
recurrence relation 

(2.3) xP k (x) = Pk+i(x) + b k+ iP k (x) + a fc P fc _i(a;) k = 1, . . . , n - 1, 

(Po{x) = l,P\(x) = x — bi) of the monic orthogonal polynomials 



(2.2) 



: <Pn 



m 



rn.fe_i 



m 



m.fc-i 



(2.4) P k {x) 



k 



n 



m r 



m 2 fc-i x n 



m k -i 



m 2 k-2 



associated with the vector (m l5 . . . , m 2n -i) [see IChiharal (119781 )]. These polynomials are orthogo 



rial with respect to every measure with first moments mi, . . . , Tn 2n -\ an d the recursion coefficients 
in (12.31) are given by 

(2.5) b k+1 = a + (b-a)((l-p 2k ^ 1 )p 2k + (l-p 2k )p 2k+1 ) ; k = 0,...,n-l 

(2.6) a k = (b-a) 2 (l -p 2k -2)p2k-i(l ~P2k-i)P2k ] k = l,...,n-l 

where we put p_ 1 = p = (note that a k > 0; k — 1, . . . , n). In the case T = [0, oo) the upper 
bound m~l does not exist, but we can still define for a point m^-i G Int .Mfc_i([0, oo)) the lower 
bound 

m^ = min < m k {^) /i G "PQO, oo)) with / tfdfi(t) = m^ for i = 1, . . . , k — 1 > , 

where m^ = (mi, . . . , m/c) T G Int .Mfc([0, oo)) if and only if m k > n\ . In this case, the analogues 
of the canonical moments are defined by the quantities 

nil — mT 

(2.7) Zl = l — Z = l,...,fc 

m ; _i - m l _ l 

(with itlq = 0) and related to the coefficients in the three-term recurrence relation (12 .3p for the 
monic orthogonal polynomials by 

(2.8) a k = z 2k ^iz 2k , 

(2.9) b k = z 2k _ 2 + z 2 k-i- 

Note that (12.71) defines a one to one mapping 

(2-10) fc: |IntA<„([0,oo))^(0oo)» 



from the interior of the moment space .M„([0, oo)) onto (R + ) n . Finally in the case T — R neither 
m^ nor m\ can be defined. Nevertheless, there exists also a one to one mapping 

J IntA4 2n _i(]R)— ►(RxR+) n - 1 xR 
I m2 n -i H> (6i, oi, . . . , a n _i, &„)-' . 

from the interior of the (2n — l)th moment space onto the space of coefficients in the three term 
recurrence relation (12.31) . which can be considered as the analogue of (I2.10p and is defined by 

(2.12) / x k P k (x)du,(x) = ai...a k k — 1, . . . ,n — 1, 

Jr 

(2.13) x h+1 P k (x)dfi(x) = a 1 ...a k (b 1 + --- + b k+1 ) k = 0,...,n-l, 

Jr 



[see for example IWalll (119481 )] . It should be mentioned that the relations for the coefficients in 
the three term recurrence relation (12. 3p could be read backwards, because any measure with 
support on an interval [a, b] or on the interval [0, oo) is also a measure on the real line. In 
other words, a measure on the real line is supported on the interval [0, oo) if and only if the 
coefficients in the recurrence relation (I2.3P for the corresponding monic orthogonal polynomials 
allow a representation of the form (I2.8P and (I2.9P with non- negative values Zk- Similarly, the 
measure is supported on a compact interval if and only if a representation of the form (12 .5p and 
( 12. 6 p holds, where the quantities p k vary in the interval [0, 1]. In the following sections we will 
use the canonical moments and corresponding quantities on the interval [0, oo) and the real line 
for the definition of distributions on the corresponding moment spaces. 

2.2 Distributions on unbounded moment spaces 



In Sec tion 4 we will show that some of the results of lChang et al.l (119931 ) and lGamboa and Lozada-Chang 



(120041 ) hold for a rather broad class of distributions on the moment space Ai n ([a, b]). For the 



definition of this class, let for k > 1 

f k : (0, 1) — > R 

be a non- negative integrable function with J Q f k (x)dx > 0, then a distribution on the interior of 
the moment space M. n ([a, b]) is defined by 



(2.14) f n (m n ) = ]~[c k , n f k (p k (m n ))l {m - 



<m,i.<m7"}' 



fe=l 



where p k (m n ) is the fe-th canonical moment defined in (12. ip and 

f k {x){x-x 2 ) n - k dx\ 

'0 



c k ,n = ({b-a)^ n+ ^ 2 f k (x)(x-x 2 ) n - k dxY (k 



n) 



are normalizing constants such that f n is a density on A4 n {[a, b}) (see the proof of the following 
Theorem 12. ip . Our first Theorem gives the distribution of the canonical moments corresponding 
to the random vector rrin with density f n defined in (I2.14p . 

Theorem 2.1. Suppose that m n is a random vector on the moment space Ai n ([a, b}) with density 
f n defined in (I2.14p . Then the canonical moments p^rrin) , . . . , p n (m n ) are independent and 
p k (mn) has the density 

c k , n (b - a)< n+1 ^ 2 f k {x){x - x 2 ) n -H {Q<x<1] 

for 1 < k < n. 



Proof: The fc-th canonical moment p k depends only on the first k moments mi, 
implies for the Jacobian determinant of the mapping ip n defined in (12. 2p 



, rrik, which 



(2.15) 



dip r 



dm n 



n 

fc=i 



dm k 



n< 

fe=i 



m. 



m, 



n< 

fe=i 



fe-i 



'n^ 1 -^)) 



-i 



i=l 



where we have used Theorem 1.4.9 and equation (1.3.6) in lDette and Studdenl (119971 ) for the last 
equality. Therefore, the Jacobian determinant in (I2.15P simplifies as 



d<fn 



dm n 



n 



fe-i 



a 



jjMi 



pi 



(6 _ a) -n(n + l)/ 2 Jj/ (1 



Pk, 



-(n-fc) 



k=l i=l k=l 

considering the product structure of /„, this gives the asserted distribution. 



□ 



In Section 4 we will show that the first k components of a random vector with density (I2.14p 
converge weak ly (after appropriat e standardization) to a normal distribution. This generalizes 
the results of IChang et al.l (Il993f ). who considered the case f k = 1 for all k > 1. For the 
construction of distributions on the unbounded moment space A1 n ([0, oo)), a special case will be 
of particular interest, that is fk(x) = x 7ft (l — x) 5k , where 7 = (7fc)fceN, $ = (^fc)fceN are sequences 
of real parameters, such that 7^, 5^ > —1 for all k > 1. In this case the density on the moment 
space JA n {[o>, b]) is given by 



(2.16) 

where 
(2.17) 



./, 



(7,<5), 



m n 



.[0,6] 



n 

k=l 



— \ 7fc 

m k -m k x 



m, 



m. 



mi 



m k 



m k ~ m k 



<h: 



1 



{m k <m fc <m+}' 



,[a,b] 



|(6-a) n(n+1)/2 /" x n - fc+7fc (l-a:) n " fc+5fc ^} ' 



i s the normalizing constant. The choice of the density ( I2.16P is motivated by results of lDette and Studden 
( 119951 ) who showed that the empirical distribution of the (appropriately normalized) roots of the 



Jacobi polynomials P k ' k (x) converges weakly to a distribution with unbounded support if 
jk — > 00 or 5k — > cxd. Note that for the density fn' the canonical moment p k has a Beta 
distribution Beta(jk + n — k + l,5 k + n — k + 1). In the following we use densities of the form 
(I2.16P to construct a distribution on the unbounded moment space 



(2.18) 



M n ([0, 00)) = {rrinifi) = (m 1 (/i),...,m n (/i)) T |/iGP([0,cx)))}. 



For this purpose, recall that the relation ( 12. 7p defines a one to one mapping between the moment 
space lntM n ([0,oo)) and (M+) n . 



Theorem 2.2. Let fn ' denote the density defined in (12. 16[) on the moment space /Vf ra ([0, d}) 
corresponding to the probability measures on the interval [0,d], where the parameter sequences 
^( d ) — (^ ^) fceN) <£(<*) = (Sj. )ken depend on length d and satisfy 

if ► lk > -l, -^ ► S k e R + 

d— >oo a d— >oa 

for all k > 1 . Then for d — > oo the density fn ' converges point-wise to the function 



(2.19) 9 r<^) = 4°-»f[( '"'""i Ye, P U "*- m i )l, 






i 

= c| l °' oo) ]^[2; fc (Tn n ) 7fc exp(-(5 fc z A; (m n ))l {2fc(mn)> o}, 
fe=i 

where the constant ci zs gzi>en fry 



c [0,oo. _ 1 f °_k_ 



n 



Moreover, g„ defines a density on the unbounded moment space .A/f n ([0, oo)). 

Proof: The fact that g4 7 ' is a density is evident from the transformation in the proof of Theorem 
12.31 below, we prove here only the convergence. For a fixed point ra^ e .M n ([0, oo)), there exists 
a4 € N with m„ e _M n ([0,<i]) for all d > do. Let p n (m„) denote the vector of canonical 
moments corresponding to the vector rrin in the moment space .M n ([0, d]). We will show at the 
end of this proof that 

(2.20) p k (m n ) = M^i(i + (i)), k = l,...,n, 

d 

where the quantities Zk(mn) are defined in (12.71) . Observing this representation and the definition 

(j2.16p . it follows for d ->■ oo 

n 

(d) , . (d) 



d -H >+-+^ >) C M J J^ (mn )7 fc eX p(-4^(mn))(l + o(l)). 



fe=i 



Finally, we obtain from (" 12 . 1 T[) for the normalizing constant 

d -h[ d) +-+7i d) ) c jM = rf -n(„+i)/2-( 7 ( d ' + ... +7 ( d ») A , r( 7 i d) + 5 { k d) +2n-2k + 2) 



„( d ) i ^ z, i 1 \r(x( d ) 



\ T( lk a > + n-k + l)T(5 k a> +n-k + l) 

d -n(n + l)/2- { ^ + - + ^) TT & ) 7 (1 + o(1) ) 

ct°°\l + o(l)), 



which proves the assertion of the Theorem. For the remaining proof of the representation (12.201) . 
let fj, be a measure on the interval [0, d] with first moments given by m n and let v denote the 
measure on the interval [0, 1] obtained from \x by the linear transformation x i— > x/d. We write 
Pk(/J>) for pk(mn) and Zk(fi) for z k (mn). Invariance of the canonical moments under linear trans- 
formations yields Pk(/J>) = Pk(y). The recursion variables of the measure v can be decomposed 
as 

(2.21) z k {y) = (\-p k ^{y))p k {y). 

The Stieltjes-transform of [i has a continued fraction expansion 

dfi(x) 1 I z 1 (jm)z 2 (jm) I z 3 (/j,)z 4 {h) I 



(2.22) 



o C-x \C-Zi(fi) C- (z 2 (fi) + z 3 (fi)) \(- (zM + z 5 (fi)) 



for (6C\[0,d]. Since the measure \x is obtained from v by linear transformation, the contin- 
ued fraction expansion of the Stieltjes-transform of \x can be written in terms of the recursion 
coefficients of v, 

dfi(x) 1 d 2 zi(v)z2(is) d 2 ''z$(y)zi(y) 



o (-x \(-dzi(v) \C,-d(z 2 {u) +z 3 (u)) \(-d(z A (v) + z 5 (u)) 



see Theorem 3.3.3 in lDette and Studdenl (J1997I )]. A continued fraction expansion as in ( 12. 22ft is 



unique, which yields dz k {y) = z k (fi). With equation (I2.2ip we obtain 

z k {v) 1 Zk{n) 



Pk(fi) = Vk{v) 



l—Pk-i(v) dl-pk-i(fj) 
for k > 1. The first canonical moment is given by 

mi-mj" mi zi(m„) 

Pi{ma) = — = -r = — -, — • 

mj — m 1 a a 

and equation (I2.20p follows by an induction argument. □ 

The following Theorem is essential for the asymptotic investigations in Section 4 and gives the 
distribution of the the vector z n = (zi, . . . , z n ) T corresponding to a random vector point on the 
moment space .M n ([0, oo)). 

Theorem 2.3. Let m n E .M n ([0, oo)) be governed by a law with density gn' , then the recursion 
variables Zn = ip n (mn) defined by (12. 7p are independent and gamma distributed, that is 

z k ~ Gamma( , y k + n — k + 1, S k ) k — 1, . . . , n. 



Proof: By its definition in (\2.7\i , the random variable Zk depends only on the moment m l5 . . . , rrik, 
therefore the Jacobi matrix of the mapping ip n is a lower triangular matrix. We obtain for the 
Jacobian determinant 



dm n 



dz n 



n ft 



n—k 

c 

fe=l fc=2 k=\ 



H 



where the third identity follows from the definition of the Z{ in ( 12. 7p . Considering the second 
representation of the density in (j2.19p . this yields the claimed distribution. □ 

We will show in Section 3 that densities of the form (" 12. 16[) arise naturally as distributions of 
moments corresponding to random spectral measures. We conclude this section with a discussion 
of distributions on the moments space corresponding to measures on R. For the sake of brevity 
we restrict ourselves to moment spaces of odd dimension, that is 

A*an-i(R) = {m 2n _ 1 (/i) = (m 1 ( ; u),...,m 2n _ 1 (/i)) T | /iGP(M)}. 

To derive a class of distributions on .M 2n _i(IR) we consider the moment space M.2n-i{[~ s, s]) 
with s — > oo and a density of the form ( I2.16P with parameters varying with s. The proof of the 
following result is similar to the proof of Theorem 12.21 and therefore omitted. 

Theorem 2.4. Denote by fj^_{ the density defined in (I2.16P on the moment space M.2n-\{{— s, s]), 
where the parameters satisfy 

Tali = S2k-is 2 + o(l), <5 2 Vi = ^_ lS 2 + o(l), 
12k = 7* + o(l), 5$ = 5 2k s 2 + o(s 2 ) 
with 7fc > — 1, 5k > 0. Then / 2 n-i converges point-wise to the function 



(2.23) h^\{m2n-i) = f[ V — ex P H2k-ib 2 k (nhn-i)) 

fc=i 7r 



«-l r -y k -\-2n~2k 

X II r(^\ 9n - 9 ^ °fc fc ( m2 "- 1 ) eXP (-^2feQfc(m2n-l)) l{a fc (m a „_i)>0}- 



fc= A T( 7fc + 2n - 2fc) 

Moreover, the function /4n-i defines a density on the moment space .M 2n _i(IR). 

The following result is the analogue of Theorem 12.31 

Theorem 2.5. Let vno,n-\ £ A^n-iO^) ^ e a random vector with density ^n-i defined in (j2.23p . 
T/ien i/ie random recursion coefficients (b±, a\ . . . , a n -i, &n) T = £271-1(^20-1) in i/ie recurrence 
relation ( 12. 3p /or £/ie orthogonal polynomials associated with m 2n _i fry ( 12. 4p are independent and 

afc ~ Gamma{~fk + 2n — 2k, 62k)- 
9 



Proof: Suppose \i is a measure with first moments given by rri2 n _i and let Pi(x), . . . , P n (x) de- 
note the corresponding monic polynomials. The recursion coefficients can be calculated by (I2.12p 
and (I2.13p . This implies that the coefficient b k depends only on the moments mi, . . . , m 2k -i, while 
a k depends only on mi, ... , m 2k . Consequently, the Jacobi matrix of (, 2n -i is a lower triangular 
matrix. From (I2.12p we obtain 



dm 2k 







da k da k 
and by the identity (12.131) it follows 

dm 2 k-i d 



x k P k {x)d[i{x) = ai . . . a k . 



db k db k 

The Jacobian determinant is therefore given by 



— - / x P k - X {x)dix(x) =a 1 ... a fc _i. 



dmzn- 



d^2n-l(m2 n -l) 



n-l 

n 

k=l 
n-l 



dm 2k -i dm 2 k \ dm 2n - 
dbk da k J db n 



n-l 



:Q]\ai...a fc _i) 2 j(ai...a fe _i) = JJ 



,2n-2fc-l 



k=2 



k=l 



X~f,S) 



Considering the product structure in the density /4n_i m terms of the recursion coefficients, the 
assertion of the theorem follows. □ 



3 Random spectral measures 

Let A be a linear, self-adjoint operator on a n-dimensional Hilbert space % with inner product 
(•, •) and cyclic vector e\ £ % (i.e., ei,Ae\, . . . ,A n ~ 1 ei are linearly independent). The spectral 
Theorem yields the existence of a unique probability measure /i on the real Borel field, such that 



(3.1) 



(e 1 ,A n e 1 ) = m n {n) 



for all n > 1 [see iDunford and Schwartzl ( 119631 )] . This measure defines the unitarily equivalent 
L 2 -space in which the operator A is represented by the multiplication f(x) h- >■ xf(x). We call /i 
spectral measure of the matrix A. If Ai, . . . , A n denote the eigenvalues of A and u%, . . . , u n is a 
corresponding system of orthonormal eigenvectors, the spectral measure can be written as 



(3.2) 



/' 



y^wA 



i=\ 



where the weights are given by lOj = |(Mj,ei)| 2 and S x denotes the Dirac measure in the point 
x. The identity (13. 2p follows easily from the fact that the matrix with columns u\,...,u n 



10 



diagonalizes A. If the eigenvalues are distinct, lie in the interior of T G {[a, b], [0, oo), R} and the 
vector of weights (wx, . . . , w n ) T is contained in the simplex 



Sim r 



,Wx, 



w 



n; T e 



n 

Wi > 0,^2w { = l| 

i=X 



then the first 2n — 1 moments of the spectral measure satisfy m>>.„.--\ (p) G Int M.2n-i(T) which 
follows by an application of Theorem 1.4.1 in iDette and Studdenl (119971 ). 



We consider in this section random spectral measures associated with three central distributions 
in random matrix theory: the Hermite (or Gaussian), the Laguerre and the Jacobi ensemble. 
These classical ensembles are distributions on the space of self-adjoint matrices A G K nxn with 
real (K = M), complex (K = C) or quaternion entries (K = H). The Hermite ensemble arises in 
physics and is the distribution of the matrix 

A=^=(X + X*), 

where all real entries in X G K nxn are independent and standardnormal distributed and X* 
denotes the adjoint of X. The eigenvalues of the Hermite ensemble have the joint density 



(3.3) 



/ G (A) = c G |A(A)|^n e " 



Af/2 



where (3 = dim^ I K, A (A ) = | X <7 -(-\> — A») is the Vandermonde determinant and cq is a normaliz- 



ing constant [see iMehtal (120041 )] . The Laguerre ensemble appears in the study of singular values 



of a Gaussian matrix, it has the eigenvalue density 

n 

(3-4) f L (X) = cl\A(X)f\l^e- x 'l {Xi> o } , 



i=X 



with parameter a > 1 and a normalizing constant c a L . Finally, th e Jaco bi ensemble is motivated 
by multivariate analysis of variance [MANOVA; see iMuirheadl ( 119821 )] and is defined by the 
eigenvalue density 



(3.5) 



/j(A) 



a.b\ 



A(A)|'nW-A*) 6 l 



{0<A,<1}- 



4 = 1 



a,fe 



with a,b > — 1 and a normalizing constant cj . It is a common feature of all three classical 
ensembles that the matrix of the orthonormal eigenvectors is Haar distributed on the group 
of ort hogonal (unitary /symplectic) matrices and independent from the eigenvalues [see iDawid 
( I1977I )]. As a consequence, the matrix distribution is uniquely determined by the eigenvalue 
density. Since the eigenvector matrix is Haar distributed, the first row of the matrix is Haar 
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distributed on the unit sphere in K n and the weights Wi, . . . ,w n follow a Dirichlet distribution 
with density 

(3-6) T(l> l ■■■ Wn 1 i^s imn }- 

The identity for the spectral measure in (13.21) motivates the following definition of distributions 
on the set of probability measures. 

Definition 3.1. Let 

H = w 1 5 Xl H h w n 5 Xn 

denote a probability measure with random support points A = (Ai, . . . , A ra ) and random weights 
with W\, . . . , w n with Dirichlet density (j3.6p . where the weights and (random) support points are 
independent. 

(i) If the density of A is given by ( I3.3J) . we call \x random spectral measure of the Gaussian 
ensemble, i.e. \i ~ G(3E n . 

(ii) If the density of A is given by ( 13. 4ft . we call \i random spectral measure of the Laguerre 
ensemble, i.e. jjl ~ Lf3E n (a). 

(Hi) If the density of A is given by ( 13. 5p . we call \x random spectral measure of the Jacobi 
ensemble, i.e. \i ~ J(3E n (a,b). 

Note that these distributions are well defined for all > 0. If /3 G {1,2,4}, we obtain the 
spectral measures of matrices from the classical ensembles. The results in this section show that 
there is a connection between classical ensembles and the distributions on moment spaces as 
given in Section 2. More precisely, the moments ( 13. ip of the spectral measures are distributed 
according to the densities defined on A1 ri ([0, 1]) and obtained on the unbounded moment spaces 
by Theorem 12.21 and Theorem 12.41 Our first Theorem considers the spectral measure of the 
Jacobi ensemble. The proof relies on a tridiagonal matrix mo del for the J(5E n {a : b) dis tribution 



for all values of /3 > 0, which was proved in a recent paper by lKillip and Nenciul (120041 ) . 



Theorem 3.2. Let \x ~ Jf3E n ('y Q ,6o) denote a random spectral measure of the Jacobi ensemble, 
then the distribution of the corresponding random moment vector m2 n -i(/-t) on M.in-\ ( [0, 1]) is 
absolute continuous with density / 2 n-i defined in ( I2.16P where the parameters of the density are 
given by 

72fc -i = (| - 2)(n - k) + 70 ; 5 2k -i = (| - 2)(n - k) + 5 
for 1 < k < n and 

72k = (f - 2)(n - k) ; S 2k = (f - 2){n - k - 1) + 7o + 6 
for 1 < k < n — 1 . 
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Proof: It follows from Theorem 2.2 and Proposition 5.3 in iKillip and Nenciul ( 12004 ) that /i is 
the spectral measure of the tridiagonal matrix 



M 



J >) 



c\ 



ci d 2 



\ 



\ 



' • c n _i 

Cji— 1 Q>n J 



where the entries are given by 

dk =2?2Jfe-2(l ~P2k-3) +P2fc-l(l ~P2k-2) 
Ck =VP2k-l(l ~P2k-2)P2k(l -P2k-l) 

with p_i = po = and pi, ■ ■ ■ ,P2n-i ar e independent random variables distributed as 



Pk 



{Beta ( 
[Beta( 



2n-k n 2n-k-2 



l P 



4 
2n-k-l 



/? + 7o + <^o + 2) , fc even, 



/3 + 7o + l,^^i/3 + 5 + l), fcodd. 



Note that IKillip and Nenciul (J2004J ) define the Jacobi ensemble by the eigenvalue density 

n 

c\A(X)f J](2 - ^) 70 ( 2 + ^) 5o M-2<x,<2} 

and therefore the matrix found in their paper is the transformed matrix 4J n — 2I n . Additionally, 
they work with a beta distribution on the interval [—1,1] which is obtained from the usual beta 
distribution on [0, 1] by the transformation x H- 1 — 2x. 

The tridiagonal matrix J n defines monic polynomials P\(x), . . . , P n {x) via a recursion ( 12. 3 p with 
recursion coefficients 

bk — dk (1 < k < n), a>k = c 1 (1 < k < n — 1). 

Indeed, the polynomial Pk{x) is the characteristic polynomial of the upper left (k x A;)-subblock of 
the m atrix J n . The orthogonality measure of these polynomials is precisely the spectral measure 
fj, [see iDeiftl (120001 ) for the corresponding statement for orthonormal polynomials] . Therefore, 
the recursion coefficients in the recursion (12. 3p of the monic polynomials orthogonal with respect 
to the measure [i are given by 



h =P2fe-2(l ~P2k-z) +P2fc-l(l ~P2k-2) ; k = 
dk =P2k-li} -P2k-2)P2ki} ~P2k-l) ; k = 1, 



l,...,n, 

. . . ,71 — 1. 
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By identity (J2.5P and (12. 6p . P2n-i = (pi, ■ ■ ■ ,P2n-i) T is exactly the vector of canonical moments 
of the spectral measure /i and by definition, their joint density is given by 



=1 

fc=l 

n 

=1 
n-l 

n^r^-'ci-p-; 



fc=l 

n-l 

(2n-2fc)/3/4-l n _ ti ^ i ^(2n-2fc-2)/3/4+ 70 +<5o + l 

fc=l 

n 

fe=l 

n-l 

(n-fc)/3/2-l/ 1 _ \(n-fc-l)/3/2+7o+5 +l 



X 

fc=l 



Since the eigenvalues of the matrix J n are contained in the interval (0, 1), the moments m^n-iin) — 
V 9 2n 1 -i(f2n-i) of the spectral measure are in the interior of the moment space A^2n-i([0, 1])- The 
Jacobian of the transformation y^n-i is given by rifc=i(Pk(l — Pk)) 2n ~ l ~ k , which gives for the 
density of the random moments rri2 n -i(fj,) 

fm(man-i) ^Hp^m^)^-^/ 2 -^^^ 
fe=i 

n-l 

xn^(^- 1 )^ W2 - (2 "- 2fc - 1) - 1 (l-^(m 2n _ 1 )) (n - /: - lW2 - (2n - 2fc - 1)+70+ ^ 1 
fc=i 

n 

=cnm-i(m 2n _ 1 )W 2 - 2 )("- fe )+^(l-P2 fe -i(m 2n _ 1 ))^ 2 - 2 )("- fe ) +5 ° 

fc=i 

x ffj*(fi^i)<^<"-*>(l -p2 fc (m2n-i)) W2 - 2)( "- fe - 1)+70+5 °. 
fc=l 

This is a density as in (12.161) with the asserted parameters. □ 



An interesting case is obtained by the choice /3 = 4. Here the parameters 7^ 5k of the density 
fzn-i do n °t depend on fc. If additionally 70 = 5 = 0, the vector of moments m^n-i is uniformly 
distributed on A^2n-i([0, 1]). In other words, starting from a moment vector m2n-i drawn 
uniformly from the moment space, we obtain a random measure 



^ = Yl Wk6x " 



fc=i 
with first moments given by Wan-i an d with support points distributed according to the density 

(3.7) /A(A) = |A(A)| 4 l {0 <A fc <iv fe} . 
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The roots of the n-th orthogonal polynomial P n (x) with respect to the measure \x are precisely 
the support points Ai, . . . , A n . Consequently the roots of the n-th random orthogonal poly- 
nomial P n (x) as s ociate d with m 2n _i. are distributed according to the density (13.71) [see also 
Birke and Dettel ( boOflf l]. 

The following results show that the spectral measures of the Laguerre ensemble and Hermite en- 
semble have moments distributed according to the density obtained in Theorem 12.21 and Theorem 
2.4[ respectively. The proofs are analogous to the proof of Theorem 13.21 and use the tridiago- 
nal m atrix models for the Laguerre and Hermite ensemble provided by iDumitriu and Edelman 
(booj ). 



Theorem 3.3. Let \x ~ Lf5E n {^o) denote a random spectral measure of the Laguerre ensemble, 
then the distribution of the corresponding random moment vector m2 n _i(u) on ,M2n-i([0, oo)) is 
absolute continuous with density gn' defined by f !2.1Q[) . where the parameters of the density are 
given by 5k = —1 for 1 < k < 2n — 1 and 



72fe-i=(f-2)(n-A;)+7o 
72fc=(i-2)(n-fc) 



1 < k < n, 
1 < k< n-1. 



Theorem 3.4. Suppose \x ~ G/3E n is a random spectral measure of the Gaussian ensemble, 
then the distribution of the corresponding random moment vector m2 n -i(/-0 on the moment space 
A^2n-i(^) i s absolute continuous with density /i 2 n-i defined by (I2.23J) where the parameters of 



the density are given by 



S 2k -i = \ (1 < k < n), 



>2k 



--1 (l<k<n-l), 



and 7fc = (| - 2)(n - k) - 1 (1 < k < n). 



4 Weak convergence of random moments 

In this section we study the probabilistic properties of random vector on the moment spaces 
M. n ([a, b]), .M„([0, oo)) and .M n (M) distributed according to the measures introduced in Section 
2. We begin with random momen ts defined on t he m oment space corresponding to probability 



measures on a compact interval. IChang et al.l (119931 ) and iGamboa and Lozada-Changi (120041 ) 



investigated the uniform distribution on M. n ([a, &]), and we first demonstrate that weak conver- 
gence of random moment vectors can be established for a rather broad class of distributions on 
M. n ([a,b}). An important role in the discussion of moment spaces corresponding to probability 
measures with bounded support [a, b] plays the arcsine distribution v with density 



dv{x) 



TV 



^/(x-a){b-x 



:1 



{a<x<b} 



dx. 
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The c anonical moments of the arcsine distribution are given by 1/2 [see iDette and Studden 
(119971 )] and therefore its sequence of moments could be considered as the center of the moment 
space JA([a, b]). The following statements establish the asymptotic properties of the (random) 
canonical moments corresponding to distributions on the moment space Ai n ([a,b]) defined in 
(j2.14p . The only assumption necessary to obtain almost sure convergence is that for all e > 

(4.1) f 2 f k {x)dx > 0. 



2- £ 



Throughout this paper the symbol — > stands for weak convergence. 

Theorem 4.1. Suppose that the distribution of the random moment vector m n e Ai n ([a,b\) is 
absolute continuous with density f n defined in ( 12.14ft , where the functions fk satisfy for all e > 
condition (14.11) . and denote by p k n the k-th canonical moment of m n (k — 1, . . . ,n). 



(a) Ifn—¥ oo ; then almost surely for any k > 1 

In) 

Pi — ■> 



(«) , I 

2' 



n— >oo 



(b) If additionally the function fk in the density (I2.14H is continous at | and fk{\) > 0, then 
the k-th canonical moment corresponding to m n satisfies 

v^^-D-^^ci). 

n— yoo 

Proof: For notational convenience, we consider only the case [a, b] = [0,1]. We introduce the 
random variable q k = p k — ~ with density 

Ck,nfk(x + 2X4 — X ) U ~ l{-l/2<a«l/2} 

and show q k n — >■ 0, which proves the first assertion of the Theorem. For < £ < |, let 
U(e) = [s,e], then clearly, 



1 = c Kn / f k (x + |)(| - x l ) n - k dx + c k , n / f k (x + |)(| - x^^x, 

JU(e) JU{e) c 

and we obtain 



( f r 1 \\,\ 2\n-kj \ 1 1 JU(e) J^ X + 2H4 

C fc , n / /fc(x+| i-X 2nfc rfx = 1+ , t 

^ (|) / fe (xH-|)(i-x 2 r- fc dx 






f u(£) cfk(x+m-x*y-*dx 






n— fe 
^ - Z 2 ) n ~ k Iu(e)C fk(x + \)dx ~ \j=&) C£J ' 
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» 



P q reU(e 



,c 



*,/' 



where c e j is a positive constant by condition (14.1 p and independent of n. This yields 

r / 1 _ r 2 \ fe ~ n 

/ / fc (x + i)(i-x 2 r^x< -- r 

which implies 

1 - £ X 



n=fe n=k ^ ' 

By the Borel-Cantelli Lemma, q k converges to almost surely, which proves the first part of 
the Theorem. 

For a proof of part (b) we note that the rescaled canonical moment y8n{p k — ~) has the density 



(4.2) 



Cfc, n ^=2 '"' /fe(| + ^=) (l - —J 1{- V ^<x<v / 2TI}' 



\/> 



and the assertion follows from the convergence of the density (14.21) to the density of the standard 
normal distribution [see the convergence Theorem by IScheffa (119471 )] . Clearly, 



2 \ n—k 



( X 2 \ 



IU-^ 2 /2 



while the factors in (14. 2 p not depending on x have the integral representation 

With the definitions 

I ■= [ f{\)e- x2/2 d X) I' n := [ ^ f k {\) (l - f) dx, 

Jr J-V2T1 \ zn J 

we have the inequality 

(4.3) \i n -i\<\i n -i' n \ + \r n -i\, 

and we need to show that the right hand side converges to 0. Obviously \I' n — 1\ converges to 
by theorem of dominated convergence. For the first term on the right hand side of ( 14. 3 j) we have 
for n sufficiently large 

/\/2n pi 

_\fk( l 2 )-fkC2 + -7t)\e- x2/2 dx= / \f k (l)-f k (l + l)\V2^e- nx2 dx. 
-J2H J-i 
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For given £>0we choose 5 > 0, such that \fk(h) ~ fk(h + f )l < e f° r & U x £ ( — $, $) then 



\In-l' n \<ef V2^e~ nx2 dx+ [ |/ fc (l)_/ fc (I + |)|v^e- 

~/(-(5,<S) J(-8,8) c 

<e I ' e- x2 l 2 dx + / |/ fe (|) - / fc (| + f )|v / 2^e~ na;2 cix 

and the second integral converges to by the theorem of monotone convergence. 



dx 



D 



The weak convergence of the canonical moments was shown by IChang et al.l ( 119931 ) for a uniform 
distribution on .A/f n ([0, 1]) (i.e., fk = 1). These authors showed weak convergence of the vector 
m^ of the first k components of a uniformly distributed moment vector rrtn 
the moment space .M n ([0, 1]), that is 



(mi, 



, m n ) on 



(4.4) 



l-l/'^W 



v 



8nA-\m^ - m k {v)) -^ N{0,I k ) 



where rrikiy) denotes the vector of the first k moments of the arcsine distribution and A is a 
k x k lower triangular matrix with entries 



i-,3 



-2i+2 



2i 
i-j 



3 < i- 



By part (6) of Theorem 14.11 it is easy to see that the weak convergence in f)4.4j) holds for the more 
general densities f n on J\4 n ([0, 1]). We conclude this paper with a discussion of corresponding 
results for distributions on the non compact moment spaces J\4 n ([0,oo)) and J\4 n (M.). In this 
case the analogs of the arcsine distribution in this context are the Marchenko-Pastur distribution 
defined by 



\/x(4: Xi 

dr)(x) = l{o<x<±}dx 

and Wigners semicircle distribution on the interval [—2, 2], that is 

1 



(4.5) 



dp(x) 



2n 



V4-X 2 ! 



{-2<x<2} 



fix- 



sec 



Nica and Speicherl (120061 )). The moments of the Marchenko-Pastur law r\ are the Catalan 

m n (r)) = c n 



numbers c n defined by 
(4.6) 



2n 



n+ 1 \ n 
and the moments of the semicircle law p are given by 



n6N, 



m n (p) 



—^r-r ( m ) if n = 2m 
if n = 2m — 1. 
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Our next results establish the asymptotic properties of the quantities z k corresponding to a 
random vector on the moment space A^„([0,oo)) with density gj?' defined in (12.191) . The 
following result is a well known consequence of the asymptotic behavior of the density of the 
Gamma distribution and the proof therefore omitted. 

Theorem 4.2. Suppose m„ is a random vector of moments on the moment space A4 n ([0,oo)) 
with density g„ ; where the 7^ are fixed, 5\ = ■ ■ ■ = S n = n, and let z£ denote the k-th 
component of the vector z, n = (z[ , ...,%)■ Then the standardized random variable z^ is 
asymptotically normal distributed, that is 

v^(4 n) -l)-^A/-(0,l). 

n— >oo 

By Theorem 14.21 the vector 

v^K (n) -i) = v^((4 n) ,...,4 n) ) T -(i 5 ---,i) T ) 

is asymptotically multivariate normal distributed. In order to derive a corresponding statement 
of the random vector m^ = ip^ l (z^ ) we will use the Delta method and study first the image 
of the vector (1, . . . , 1) T under the mapping ip^ 1 . 

Lemma 4.3. Let (c ra ) ra >i denote the sequence of Catalan numbers defined in (14. 6p . then 

V>n(Cl,...,C n ) = (1,...,1) T . 

Proof: The proof presented here relies on the combinato rical interpretation of the Catalan 



numbers and a recursive algorithm given in ISkibinskyl (1l968h to calculate the moments in terms 
of the variables z k . The k-th Catalan number counts the paths in N x N starting in (0,0) and 
ending in (2/c, ), where one is only allowed to make steps in the direction (1,1) or (1,-1). 



Skibinskyl ( 119681 ) defines the triangular array {gi,j}ij>o by g it j = for i > j, g j = 1 and the 



recursion 

(4.7) g itj = g itj _i + Zj^+ig^u, 1 < i < j. 

He showed that g kj k = m^. Consequently, if zi — 1 (i = 1, 2, . . . ) the quantity g k ,k is the number 
of paths through the lattice {(i, j)}ij>o, starting in (k, k) and ending in (0,0), where in each 
vertex we can only make steps upward or to the left and where we are not allowed to cross the 
diagonal {(i,i)}. This number is exactly the k-th Catalan number c&. □ 



Theorem 4.4. If the vector of random moments m„ G .M n ([0, 00)) is governed by a law with 
density gn' , where 8± = ■ ■ ■ = S n = n and the 7^ are fixed, then the projection rri^ 1 = Il£(m„) 
of rrin onto the first k coordinates satisfies 



V^C~\mP - m k ( v )) -^ A4(0, 1 



k) 



19 



where the vector m k (rj) = (c\, . . . ,c k ) T contains the first k moments of the Marchenko-Pastur 
distribution and C is a lower triangular matrix with entries C\ y \ = ■ ■ ■ = c k)k — 1 ; an d 

^ = C-%)-G-?-i) i<L 

Proof: It suffices to calculate the Jacobi matrix 



dz k 



k. 



of the mapping ip k 1 at z® = (1, . . . , 1) T , then the independence of the recursion variables z k and 
Theorem 14.21 yield the convergence 

v^(mf n) - m k ( V )) = ^C{z[ n) - z£) + o P {\) -?-> Af k (0, CC T ). 

n— >oo 

Note that the moment vfi{ depends only on z\, . . . , zi and consequently C is a lower triangular 
matrix. To identify the entries of the matrix C we consider the triangular array {gi,j}i,j>o defined 
in (14. 7ft . For a fixed r with 1 < r < k we introduce the notation 

and obtain a new triangular array {ui,j}i,j>o- Obviously we have Uij = for i > j and the other 
values of Uij are determined by the initial condition it j = and the recursion 

(4.8) u^ = Uij-! + itj-ij + 8 rjj -i +1 g?_ ld , l<i<j, 

were 8ij denotes the Kronecker symbol and g>°- is the coefficient in the recursion (14.81) . if all Zi 
are equal to 1, that is 

o fi + J\ ( i + J 



The numbers gfj are sometimes called generalized Catalan numbers [see iFinucanl (1l976l )]. At the 



end of this proof we will show that the entries in the new triangular array are given by 

With this identity we obtain for the entries of the matrix C 

(zS) = Ui 



oz r \i — rj \i — r — I 

for 1 < r < i, which proves the assertion of the Theorem. 
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The proof of (14.91) follows by an induction argument for the row number i. Obviously we have 
u j = and therefore U\j = as long as j < r. The value of u\, r is equal to g$ r = 1 and for 
j > r we have u\j = 1. These values are also given by the formula (14. 9p . In the induction step 
i — 1 — > i we perform an induction for j and have to distinguish between five cases with different 
recursion for Uij. The scheme in (I4.10p illustrates these cases in the i-th row. 

(4.10) ' • . (i, i) ' ■ . (i, i + 1) . . . (i, i + r - 2) ' • . (i, i + r - 1) ' • • (i, i + r) ' ■ ■ (i, i + r + 1) . . . 

(1) (2) (3) (4) (5) 

(1) j = i: First consider r = 1, Then case (1) is identical to case (3) and 

i (2i— 1\ /2i— 1\ (2i— 1\ /2i—l\ 

_( 2i \ _ (2i—l\ _ (2i— 1\ _ / 2i \ _ / 2t \ 
~~ U-lJ W-2/ U-3/ "" li-lJ li-2/' 

For r > 1 it is 

Ull = u^u = ( 2i - 1 ) - ( 2l -\) = ( 2l ) - ( 2l -\) - ( 2l -\) = ( 2l )-l 21 ,)■ 

l ' 1 ' x >' \ i—r J \i—r—2/ \i—rJ Vi— ) — 1/ \i— r— 2/ \i—rj \i— r— 1/ 

(2) z < j < i + r - 1: 

M i,i = M M'-1 + U *-l,i = (j-l-r) ~~ U-r-l) + ( j-r J ~~ (i-r-2) = (j-J ~~ li-r-J 

(3) j = i + r — 1: The case r = 1 was considered in (1), for r > 1 we have 



u^ = Uij-t + u^xj + g t _ ltj 

\j-l-rj \i-r-lj ~ \ i-2 ) \i-r-2J ~ \ i-l ) \ j-2 ) 

_ C»+J-1N _ ( i+i \ , fi+i-lN _ (i+j\ _ ( i+i \ 
\j—l—rj \i— r— 1/ V j— r I \j—rj \i— r— 1/ ' 

(4) j = i + r: 

«* ^-i+^-r, = (*£) - Ctl) + CX 1 ) - (E: 2 ) = (Kir) + OS) - G::ix) 

— ^ i+ J ^ _ ( i+ J ^ == < i+ i\ - ( i+ i \ 

\j-r-l) \i-r-l) \i-l) \i-r-l) 

(5) j > i + r: 

u^ = u^ 1+ u^j = e+17 1 ) - £-;) + (^ 2 1 ) - (SiJ) = £0 - G-^) 



n 



We conclude the discussion of random moments considering the moment space M.2n-xi 
Recall the bijective mapping (12.111) from the interior of the moment space A^n-iO^) onto 
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(R x R + )"' x R corresponding to the range for coefficients in the recursive relation of the 
orthogonal polynomials (12. 3p . The following results give the weak asymptotics of random recur- 
sion coefficients and moments and correspond to Theorem 14. 21 and 14.41 The proof of Theorem 14.61 
follows by similar arguments as presented in the proof of Theorem 14.41 and is therefore omitted. 

Theorem 4.5. Let the random vector m 2n _ 1 G A^2n-i(R) be governed by a law with density 
^2n-i w here 7& > —1 is fixed and 5 k = n (k = 1, . . . ,n). For fixed k denote by b k and by 
a k the (2k — l)-th component of the vector £,2n-i( m 2n-i) an d the 2k-th component, respectively. 
Then 



2nb k n) ~Af(0,l), 
2n(af-l)-^Af(0,l). 

n— >oo 

Theorem 4.6. Let the vector of random moments m2 n -i £ A^2n-i(R) be governed by a law 
with density /i 2 n-i where 5 k = n (k = 1, . . . ,n). For k e N denote by mj^ = I[ k (m2 n ^i) the 
projection onto the first k coordinates and by m k (p) = Il fe (0, c\, 0, c 2 , . . . ) the vector of the first 
k moments of the semicircle law defined in f )4.5p . Then 



2nD-\m k n) - m k (p)) -^ Af k (0, I k ), 

n— >oo 

where D is a k x k lower triangular matrix with d^j = ifi + j is odd and the remaining entries 
are given by 



di,j \ j-j I 1 i-j -I 

V 2 / V 2 

By the results in Section 3, the moment density of the three classical ensembles is the moment 
density investigated asymptotically in this Section. Although for the random matrix ensembles 
the parameters 7&, 5k depend on n, only minor changes are necessary to obtain a weak convergence 
result for the first k moments. In the case of the Jacobi ensemble, the canonical moment p k 
follows a Beta distribution with parameters behaving like ~n and therefore, 

V^(pi n) - \) -^ A/"(0, 1) 

n— »oo 

and we obtain easily the right scaling for the ordinary moments in the following Corollary. 

Corollary 4.7. Let ji n ~ Jf3E n ( r y ,5o) be a spectral measure of the Jacobi ensemble, then the 
first k moments m k (fi n ) of \i n satisfy 



yflfinA 1 {m k {n n ) - m k {v)) ► W fc (0, I k ), 

n— >oo 

where m k (v) is the moment vector of the arcsine measure and A is the k x k matrix in (14. 4p . 
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In particular, the moment convergence implies the weak convergence of the spectral measure to 
the arcsine measure. This is also a consequence of the well-known convergence of the emprical 
eigenvalue distribution to the arcsine measure, since the (unsealed) moments of the spectral mea- 
sure have the same asymptotic behaviour as the moments of the empirical eigenvalue distribution. 
By Corollary 14. 7\ the fluctuations around this limit in terms of the moments are Gaussian. A 
corresponding result holds for the Laguerre ensemble with rescaled eigenvalue density 



/ L (A)=cl°|A(A)|' 3 n^ e"^ /2 l { A l >o } . 



For the first k moments of a spectral measure fi n with this eigenvalue density we obtain from 
Theorem 13.31 together with Theorem 12.31 the asymptotic 



^/2C- l (m k (v n ) ~ m k ( V )) -^A4(0,/, 



v 

ra— »oo 



From Theorem 13.41 and Theorem 12.51 we can deduce for a spectral measure /i n of the rescaled 
Gaussian Ensemble with eigenvalue density 



/ G (A) = c G |A(A)|^n e 



n 

2, 



-PnXf/i 



i=l 



the weak convergence of the first k moments 



V/3n/2D 1 (m k {lJLn) ~ m k {p)) >N k (0,I k ). 

n— >oo 

The different scaling in these two cases comes from the scaling necessary to obtain weak con- 
vergence of the recursion parameters z k and the recursion coefficients a k ,b k , respectively. Again, 
this convergence results can be related to the convergence of the empirical eigenvalue density to 
the Marchenko-Pastur law and the semicircle distribution and give the fluctuations around the 
limit distribution. 
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